Projecting the impact of climate change on honey bee plant habitat distribution in Northern Ethiopia

Climate change significantly affects the diversity, growth, and survival of indigenous plant species thereby influencing the nutrition, health and productivity of honey bees (Apis mellifera). Hypoestes forskaolii (Vahl) is one of the major honey bee plant species in Ethiopia’s Tigray region. It is rich in pollen and nectar that typically provides white honey, which fetches a premium price in both local and inter-national markets. Despite its socio-economic and apicultural significance, the distribution of H. forskaolii has been declining, raising concerns regarding its conservation efforts. However, there is limited knowledge on how environmental and climatic factors affect its current distribution and response to future climate change. The study investigates the current and projected (the 2030s, 2050s, 2070s, and 2090s) habitat distributions of H. forskaolii under three future climate change scenarios (ssp126, ssp245, and ssp585) using the Maximum Entropy Model (MaxEnt). The results show that land use (50.1%), agro-ecology (28%), precipitation during the Driest Quarter (11.2%) and soil texture (6.1%) predominantly influence the distribution of H. forskaolii, collectively explaining 95.4% of the model's predictive power. Habitats rich in evergreen trees and mosaic herbaceous with good vegetation cover are identified as the most suitable for H. forskaolii. The spatial distribution of H. forskaolii is concentrated in the highlands and mid-highlands of the eastern and southern parts of Tigray, characterized by a colder temperature. Across the three climate change scenarios, the size of suitable habitat for H. forskaolii is projected to decrease over the four time periods studied. Predictions under the ssp585 scenario reveal alarming results, indicating a substantial decrease in the suitable habitat for H. forskaolii from 4.26% in the 2030s to 19.09% in the 2090s. Therefore, given the challenges posed by climate change, research efforts should focus on identifying and evaluating new technologies that can help the H. forskaolii species in adapting and mitigating the effects of climate change.


Study area
Tigray is the north most regional state in Ethiopia (Fig. 1), sharing borders with Sudan, Eritrea, and Amhara and Afar 29,30 .It covers an area of 54,593 km 2 with an elevation ranging from 500 to 4000 m.a.s.l.Tigray exhibits agroecological diversity, characterized by variations in temperature, rainfall, topography, soil characteristics, vegetation cover and other natural resources 31,32 .The region is characterized by lowlands (< 1500 m.a.s.l.), which account for 53% of the total area.Midland (1500-2500 m.a.s.l.) and highland areas (> 2300 m.a.s.l.) account for 39% and 8% of the region, respectively 33 .The lowlands are characterized by high temperatures, while the highlands are by low temperatures 31 .The region has an average temperature of 18 °C30,34 .Rainfall in the region is highly seasonal, with the majority of precipitation falling between June and September 32,35 , received within 50-60 effective rainy days 29 .The average rainfall ranges from about 200 mm in the northeast to over 1000 mm in the south western highlands 30 .The land use types of Tigray can be broadly categorized into bare land, agricultural land, forested land, bushland, shrubland, grassland, and arable land 36 .Similarly, major soil types in Tigray include cambisols, lithosols, acrisols, fluvisols, luvisols, regosols, and vertisols, which are characterized by a range of soil textures from light to heavy 29 .
Due to its natural landscape and diversified agro-ecology, Tigray is home to a rich biodiversity 37 .H. forskaolii (Fig. 1) is one of the most important bee plant species widely found in the region 5,18,21 .This plant belongs to the family Acanthaceae.It is a drought-tolerant perennial herb reaching a height of up to 1 m.It produces tiny white flowers with abundant pollen and nectar, its peak flowering period extends from September to November 5,18 .

Species occurrence data
Species occurrence data was collected from June to September 2020 through field surveys and discussions with key informants including bee technicians, bee experts, and beekeepers (3-5) in each village (locally Tabia) and district of the region.Once the presence of H. forskaolii was confirmed in a particular Tabia, one occurrence record with GPS coordinates was documented.To ensure comprehensive representation, each sampling point was expanded by ten additional points, following the methodology outlined by Gebru et al. 38 .To eliminate duplicated occurrence points and minimize the potential influence of environmental variations on model accuracy, we used the species distribution model (SDM) Toolbox, a Python based ArcGIS toolbox.The data were filtered to retain only one data point per grid, each 1 km by 1 km dimensions.This filtering process resulted in a total of 317 occurrence records for H. forskaolii (Fig. 2).

Environmental and bioclimatic variables
To predict the distribution of H. forskaolii across different areas, various factors including bioclimatic conditions (temperature and precipitation), soil conditions, topography and land use types were considered (Table 1).Data related to elevation and bioclimatic variables (bio1-bio19) that represent the near-current climatic conditions (1970-2000) were obtained from the WorldClim database (http:// www.world clim.org/) at a high spatial resolution, 30 arc seconds (≈ 1 km 2 ).To forecast future climatically suitable areas for H. forskaolii, projections from the HadGEM3-GC31-LL model were chosen due to its optimal performance in Ethiopian environments 39.Here, three scenarios of Shared Socio-economic Pathways (SSP) (ssp 126, ssp 245 and ssp 585) were selected across four   Additionally, soil conditions and land use land cover data were obtained from the ISRIC World Soil Information (ISRIC Soil Data Hub | ISRIC) and Land Cover CCI-LC (Geoportail UCL-Geomatics), respectively.Given that agro-ecology also plays a crucial role in a species' survival and habitat distribution, it has been included as an input variable in the model.There are several agro-ecological zonation systems that categorize the landscapes of Ethiopia 31,[40][41][42][43] .Among these, the traditional method of climate classification which primarily relies on elevation and temperature 31,43 has been widely used by policymakers and practitioners when formulating development strategies, such as agricultural land-use and natural resource management in Ethiopia 40 .There is also a strong correlation between altitude and agro-climatic zones 31 .Therefore, the traditional agro-ecological classification system has been employed in this study 41 to allow policymakers and practitioners to easily understand and apply the research output.The incorporation of these diverse environmental variables ensures a comprehensive and holistic approach to climate projections.

Variable selection procedure
Various techniques, including Spearman correlation, Principal Component Analysis (PCA) and the 'MaxentVari-ableSelection' , have been used for environmental variable selection.To reduce multicollinearity among predictor variables and prevent overfitting in the model, pairwise correlation analyses were conducted in R, employing the corrplot package.Highly correlated variables, with a Spearman correlation coefficient |R|≥ 0.8, were subsequently removed (Fig. 3).To further refine our variable selection, PCA was performed to evaluate the contribution and relationship of each variable within their respective groups of climatic, soil, topography and land use using the Stats R package.The variable with the highest contribution was chosen based on its eigenvalues and eigenvectors (Table 2) 44 .Finally, we simultaneously evaluated all variables using the R package for Maxent Variable Selection (MVS) 38 .These selection techniques helped us to identify 19 environmental variables from the initial pool of 30 variables.Then, the exclusion of variables with a zero-contribution rate (bio10 and bio11) to the model's   www.nature.com/scientificreports/prediction was performed to improve its performance.Finally, seven variables (bio4, bio8, bio9, bio17, aez, land use, texture) were selected for predicting the distribution of H. forskaolii using MaExent (Table 3).

Selection of model parameters
The prediction of the MaxEnt model can be influenced by model parameters such as feature class (FC) and betamultiplier (BM) 38 .To enhance the accuracy of species distribution predictions, adjustments were made to the feature combination and regularization multiplier using the Enmeval package.The selection of the best FC and BM combinations was based on the least Akaike Information Criterion (AICc) values [45][46][47] using the best features (linear) and regulation multiplier (RM = 3), resulting in 19 combinations.In addition, the MaxEnt iteration was set to 5000, allowing for a more comprehensive analysis.To address sampling bias, a bias file was also included in the model run.

Model calibration and evaluation
The model's accuracy is significantly improved when the training dataset includes a broader geographic area rather than relying solely on specific data 48 .To enhance accuracy and address uncertainty, we employed the tenfold cross-validation method, running it 10 times to minimize and average the results.For model calibration and evaluation, the species occurrence data were randomly split into a training sample (75% of total occurrences) and a test sample (25% of total occurrences) 49 .In this study, despite the potential loss of statistical independence between training and test data, the "bootstrap" approach was employed as a sampling technique 48,50 .The bootstrap method enabled MaxEnt to replace the entire population data with the data used for model development.The model's performance was evaluated using Area Under the Curve (AUC) and the True Skill Statistic (TSS) metrics.AUC measures how effectively a model distinguishes between presences and absences, with its value ranging from 0 to 1.A value of 1 indicates a perfect fit of the model.TSS evaluates a model's ability to detect true presence (sensitivity) and true absence (specificity), calculated as sensitivity + specificity − 1 (Eq.1).The TSS value, ranging from − 1 to 1, signifies the degree of agreement between predicted and actual points, with a score of 1 indicating perfect agreement 51,52 .

Habitat suitability classification
MaxEnt uses a habitat suitability index (HSI) ranging from 0 to 1 to evaluate the habitat's suitability for a species, considering diverse environmental variables.A value of 1 represents the highest suitability for the habitat 53 .Here, the generated habitat suitability map was categorized into four groups: not suitable, less suitable, moderately suitable, and highly habitats using the Natural Breaks method in ArcGIS and the area within each category was calculated (Fig. 11; Table 4).

MaxEnt model performance evaluation
The average area under the curve (AUC) value, based on ten repetitions of the model, is 0.906 (Fig. 4) indicating excellent predictive ability for identifying potentially suitable habitats for H. forskaolii.Additionally, the True Skill Statistic (TSS) value is 0.634, demonstrating the model's good discriminating power between points with and without the presence of H. forskaolii.Overall, the results suggested that the final model developed for predicting suitable habitats for H. forskaolii in Tigray is highly accurate, indicating its effectiveness as observed in the model performance evaluation.

Contribution of bioclimatic variables to habitat suitability
The relative importance of predictor variables and their contributions were evaluated using the Jackknife method (Table 3; Fig. 5).The land use, agro-ecology, and precipitation of the Driest Quarter (mm) (bio17) are identified as crucial factors influencing the distributions of H. forskaolii (Table 3), collectively contributing to 89.3% of the model's predictive power (Table 3).Furthermore, agro-ecology and land use have shown the highest regularized training gains (Fig. 5).These results strongly indicate that land use and agro-ecology are the most important determinants influencing the occurrence of H. forskaolii.Conversely, the contribution of environmental factors such as the mean temperature of the Driest Quarter (bio9), temperature seasonality (bio4), and mean temperature of the Wettest Quarter (bio8) was small, accounting for 4.6% (Table 3).

Environmental variable response curves
Understanding the climatic thresholds for the core habitats of H. forskaolii is crucial for effective land management and natural resource conservation efforts.Identification of these thresholds has provided insights into the suitability of the specific environmental conditions for the growth and distribution of the species.According to the response curve for land use, the natural suitable habitat for H. forskaolii mainly includes areas dominated by broad leaved evergreen trees followed by mosic herbaceous cover (Figs. 6, 7 and 8a).Whereas, areas covered with shrubs and deciduous trees that are open, covering between 15 to 40% of the land areas are less suitable habitat for the distribution of H. forskaolii.The response curve for the precipitation of the Driest Quarter (bio17) indicates that the likelihood of H. forskaolii occurrence increases with precipitation during the driest quarter and reaches its peak when the area receives about 130 mm.This suggests the importance of a certain amount www.nature.com/scientificreports/ of precipitation during the driest season for the survival of the plant.The plant species is commonly found in core areas with specific soil texture classes such as loam, clay loam and sandy clay loam soils (Fig. 8b).This is supported by the response curve for soil texture, which showed that the occurrence likelihood of H. forskaolii peaks when the soil type is loamy sand, followed by sandy loam (Fig. 8b).Additionally, the response curve for the mean temperature of the wettest quarter (bio8) and mean temperature of the driest quarter (bio9) indicated that the occurrence likelihood of H. forskaolii reaches its peak when the temperature falls below 11 °C and 8 °C, respectively.The study also revealed the highest habitat suitability in the highland and mid-highland agroecological zones compared to the lowlands.

Current habitat distribution of H. forskaolii
Under current climatic conditions, about 55.5% (30,272 km 2 ) of Tigray's land mass is identified as climatically suitable habitat for H. forskaolii (Fig. 9).These suitable habitats are mainly concentrated in the highlands and mid-highlands of the eastern and southern zones.Among the suitable areas, a significant portion was classified as highly (10%, 5467 km 2 ) and moderately (24.6%, 13,407 km 2 ) suitable habitats for H. forskaolii.These suitable habitats provide optimal conditions for the growth and distribution of the species.Most of the eastern and southern zones of Tigray are classified as moderately and highly suitable habitats (Fig. 7c).Conversely, areas less suitable for H. forskaolii, are primarily located in the central and northwest zones, covering 20.9% (11,398 km 2 ) of the region.On the other hand, large parts of western and northwestern Tigray, characterized by lowlands (Fig. 7c), are dominantly not suitable for the growth and distribution of H. forskaolii.Identifying the suitable and unsuitable habitats for H. forskaolii has important implications for the management of the species including    its conservation and cultivation.The concentration of moderately and highly suitable areas in the eastern and southern zones suggests that these areas should be prioritized for the species conservation and the white specialty honey production.In addition, honey bee population density has followed the same spatial pattern (Fig. 7a) as that of H. forskaolii distribution in Tigray.While H. forskaolii likely contributes significantly to this overlap, it may not entirely account for it since other species of bee forage plants in the area could also contribute to the high honey bee population density and suitability areas for beekeeping 17,18 .As H. forskaolii thrives in these areas, it is crucial to protect and preserve its natural habitats to ensure the long-term survival of the species and enhance white honey production in the region.

Future habitat distributions of H. forskaolii
Table 4 depicts future H. forskaolii 's species distribution predicted under three future climate change scenarios (ssp126, ssp245 and ssp585).Under all these scenarios, the predicted suitable habitat for H. forskaolii is expected to decrease across the four time periods included in this study: 2030s, 2050s, 2070s, and 2090 (Figs. 10 and 11).However, the suitable habitat for H. forskaolii is expected to decrease significantly under the ssp585 scenario and decreases exponentially with time.Under this scenario, the suitable habitat for H. forskaolii will decrease by 4.26%, 8.67%, 13.56%, and 19.09% compared to the current distribution in the 2030s, 2050s, 2070s and 2090s, respectively (Fig. 10).
In addition, a significant reduction in the suitable habitat for H. forskaolii is expected towards the end of the century under the ssp245 scenario.The result revealed that suitable habitat for H. forskaolii under the ssp245 scenario is predicted to reduce by 9.01% and 10.47% by the 2070s and 2090s, respectively (Fig. 11).Similarly, under the ssp126 scenario, the suitable habitat for H. forskaolii is expected to decrease progressively from 3.36% in the 2030s to 6.71% in the 2090s.
Under ssp585 and ssp245 scenarios areas regarded as not suitable and less suitable habitats for H. forskaolii are expected to increase, particularly in the northwestern and central zones of Tigray (Figs. 10, 11).Conversely, areas considered moderately and highly suitable habitats for H. forskaolii are expected to decrease, with the most significant reduction anticipated under ssp585 scenarios towards the end of the century.Consequently, the areas that are highly suitable habitats for H. forskaolii will be limited to a few spots in the eastern part of the region (Fig. 11).

Discussion
Climate change significantly affects the distribution, growth, and survival of many indigenous plant species, which are primary food sources (pollen and nectar) for honey bees [26][27][28] .Therefore, understanding how these plants respond to global climate change is crucial for understanding and managing their habitat distribution.H. forskaolii is one of the most important honey bee plants in Ethiopia's Tigray region, producing abundant pollen and nectar for honey bees 5 .However, due to the decline in its population distribution, further investigation into the contributing environmental factors is necessary to develop possible conservation and restoration strategies.Thus, predicting the current and future distribution of H. forskaolii under various climate change scenarios is essential for informed policymaking.

H. forskaolii species habitat distribution using MaxEnt model
The model's predictive performance was evaluated using AUC, with values ranging from 0 to 1.A higher AUC value indicates better prediction performance 54 .More specifically, an AUC value greater than 0.9 indicates highly accurate prediction 55 .In this study, an average AUC value of 0.906 was obtained indicating that the simulation www.nature.com/scientificreports/effect was excellent and, thus, the model can be used to predict the potential distribution of H. forskaolii.The result also reflected a high level of discrimination between suitable and unsuitable habitats 56,57 .The TSS values range from 0 to 1, with values close to 1 indicating higher model accuracy 52 .The TSS values were 0.634, suggesting a high probability of accurate prediction 54 , indicating the model's effectiveness in discriminating between the presence and absence of H. forskaolii species 52,54 .Overall, the results indicated that the final model developed for predicting suitable habitats for H. forskaolii in Tigray is highly reliable and accurate.

Determinant factors for H. forskaolii distribution
Under the current climatic conditions, 55.5% of Tigray's land area (30,272 km 2 ) is a suitable habitat for the growth and distribution of H. forskaolii.Moderately and highly suitable areas for H. forskaolii are concentrated in the eastern, central, south, and southeast zones, which are characterized by highland and mid-highland agroecological zones (Fig. 7c).The prevalence of suitable habitats in these zones aligns with the abundance of H. forskaolii.Notably, these areas are well known for their white honey production, primarily produced from this plant species 15 .Conversely, suitability is lower in the hot lowland agroecological zones in the northwest and west zones (Fig. 9) 31 .These results suggest that the most favorable environmental conditions for the plant's distribution are found in the highland to mid-highland areas of Tigray.Moreover, the occurrence of H. forskaolii is found to be significantly influenced by agro-ecology (Table 3), where habitat suitability is highest in the highland and mid-highland areas.These areas are characterized by lower temperatures (Fig. 7) as compared to the lowland areas in Tigray 31 .The response curve for the mean temperature of the wettest quarter (bio8) and driest quarter (bio9) aligns with these results.These results indicate www.nature.com/scientificreports/ that the likelihood of H. forskaolii incidence reaches its peak when the temperature falls below 11 °C for bio8 and 8 °C for bio9.This means that colder environments are more conducive, demonstrating the significance of temperatures in influencing the distribution of H. forskaolii.The prevalence of suitable habitats in the highland and mid-highland agroecological zones (Fig. 9), along with the crucial role of agro-ecology in predicting the plant's distribution (Table 3), suggests a potential association between preferred habitat, agro-ecology, and altitude.Thus, agro-ecology, a specific geographical area with certain environmental conditions, such as climate, soil, topography, and hydrology 58 has been a key factor in determining habitat suitability for H. forskaolii.However, the specific relationship between climate and ecological variables may vary depending on the species and environmental conditions 59 .The overall findings underscore the agroecological preferences of H. forskaolii, highlighting the necessity to prioritize specific ecological-regions in conservation and production initiatives.These suitable habitats are proposed for the conservation of H. forskaolii and could also be used as a source of planting materials for the less suitable areas in the region augmented with improved management.
Precipitation during the driest quarter (bio17) is a bioclimatic variable that indicates the total precipitation received during that specific quarter.It is commonly used to assess water availability and drought stress in a given area 60 .In our study, bio17 ranks among the top three environmental variables influencing the model's prediction (11.2%), signifying its relevance in determining the occurrence of the plant species.The response curve for bio17 indicated that the likelihood of H. forskaolii occurrence peaks between 0 and 130 mm in precipitation (Fig. 6), highlighting the importance of water availability during the driest season for the species' survival, growth and distribution.Lower bio17 values (low precipitation during the driest quarter) may also suggest that the plant species has a higher tolerance to drought stress or water scarcity.The annual precipitation in the eastern and southeastern parts of Tigray, where H. forskaolii is best suited, is lower compared to the western part (Fig. 7b).These areas also receive a small amount of precipitation during the dry season 32 .Hence, the variation in suitable habitat coverage for H. forskaolii across different administrative zones in Tigray may also be attributed to variations in rainfall distribution.
Soil texture is the most important basic physical property of soil that determines water holding and infiltration capacity 61 which in turn plays a significant role in shaping natural vegetation patterns 62,63 .In this regard, woody plant dominance is associated with coarse-textured soils, while fine-textured soil supports herbaceous species 63 .In our study, the soil textural class emerged as a crucial determining factor (6.1%) for the occurrence of H. forskaolii (Table 4).Loam, clay loam and sandy clay loam soils are suitable for H. forskaolii (Fig. 7b).While clay soil, which retain more water, is less suitable, especially in logged areas.Most of the plant's distribution is also found in clay loam and sandy clay loam soils (Fig. 8b).This indicates that the plant is adapted to low soil moisture content 64 .In semi-arid environments, coarse-textured soils like loam, clay loam and sandy clay allow water to infiltrate better during rainy season and plants with extensive root systems, such as H. forskaolii, can access this stored water during dry season 63 .
Land use plays a crucial role in shaping species distribution by directly impacting their habitats and ecosystems 65 .Changes in land use due to human activities, such as deforestation, industrial development, urbanization, and agricultural land expansion have led to declines in species abundance, diversity and ecosystem health worldwide 66,67 .The present study also identified land use as the most determinant factor influencing the distribution of H. forskaolii, contributing more than 50% to the model's prediction (Table 3).In our study, land uses dominated by evergreen trees and mosaic herbaceous with good vegetation cover that create relatively cooler microclimate have been identified as the most suitable habitat for H. forskaolii (Fig. 6).While areas dominated by shrubland and deciduous trees that shed their leaves during the hot dry season are considered less suitable for this species.Therefore, H. forskaolii is a shade-tolerant plant, well-adapted to low light conditions beneath the higher canopy of large and tall trees.Consistent with this, the response curve for the Mean Temperature of the Driest Quarter (Fig. 9), shows that the suitability decreases significantly as the temperature of the quarter increases from 10 to 25 °C.These results suggest that areas with good vegetation cover provide continues shade and a relatively colder microclimates are the natural habitat for H. forskaolii 68 .In addition, H. forskaolii exhibits adaptation to a diverse range of land use types, including riverine forests, rocky slopes, open woodlands, mixed woodlands, and wooded grasslands 5 .On the other hand, grasslands and agricultural land are made up of open spaces with annual crops and shorter vegetation at which the H. forskaolii has not developed resistance to regular grazing and trampling effects of livestock.
However, various natural and anthropogenic activities, such as frequent drought, deforestation, land degradation, agricultural land expansion, armed conflict, and other human activities, have significantly contributed to the decline in forest cover in Tigray 36,69 .Only from 1965 to 1994, the forest land coverage in Tigray decreased by 2.9% in favor of arable land and rangeland 36 .This suggests a transformation in land use patterns over those years, highlighting a trend where forest and bush areas have been reduced to make way for increased arable cultivation and grazing lands., Between 2000 and 2020, vegetation cover in Tigray increased by 2% due to a continuous reforestation and natural resource rehabilitation initiative [70][71][72] .However, during the period 2020-2022, the vegetation cover decreased by 5% of the area 72 , and forest cover declined on average by 14.5% due to the increased demand for firewood and charcoal during the wartime 69 .H. forskaolii's coverage and population growth could also be affected by military actions such as destroying trees, soil and water conservation efforts, disrupting soil health 69,72,73 .These suggest that the loss of H. forskaolii habitat is not solely an effect of climate changes but also due to other anthropogenic activities.Therefore, to revitalize and conserve the important honey bee plant, H. forskaolii appropriate land use policies, plantation and conservation practices should be implemented.Climate change can have a significant impact on honey bees by disrupting the diversity and availability of plants crucial for collecting pollen and nectar.The alteration in climate patterns poses a threat to honey bees as it destroys their natural food sources, potentially leading to adverse effects on nutrition, overall health, and productivity 28 .The current study indicates a decline in the suitable habitat for H. forskaolii, a major honey bee plant in Tigray due to climate change, particularly under the ssp585 climate scenario (Fig. 11 and Table 4).Projections under this scenario show a gradual decline in suitable habitats, from 4.26% in the 2030s to 19.09% in the 2090s (Fig. 11 and Table 4), compared to the current distribution (Fig. 10).Similarly, the suitable habitat for H. forskaolii is anticipated to decrease under the ssp245 scenario.These results underscore the negative impacts of climate change on the availability of this plant species.This could exacerbate the existing shortage of honey bee forages in the region 74 .This will also significantly affect the livelihoods of the local beekeepers by reducing white honey production, a major source of income.The scarcity of bee forage may also lead to absconding, wherein bees leave their hives in search of food, ultimately reducing honey production.Therefore, the reduction in highly suitable habitats in the eastern and south-eastern, and the reduction in moderately suitable habitats in the western, southern and central zones highlights the necessity for urgent interventions aimed at preserving, sustaining, restoring and promoting the honey bee plant populations.These areas are homes for densely populated honey bees in Tigray (Fig. 7a).Acknowledging the diverse impacts of climate change across Tigray and identifying moderately and highly suitable areas highlights the importance of prioritizing specific localities for conservation efforts aimed at preserving H. forskaolii that would be used as a source of planting material in promoting this important species through an augmented management in wider areas.

Conclusion
In this study, the MaxEnt model was used to identify potential habitats for H. forskaolii in Tigray, northern Ethiopia.Presently, about 55.5% of Tigray's total land mass is identified as suitable habitat for H. forskaolii, with a predominant concentration in the highlands and mid-highlands of the eastern, central, southeastern and southern zones.The key determinant factors influencing the distribution of this plant species include land use patterns, agroecological conditions, precipitation levels during the Driest Quarter (mm) (bio17) and soil texture, collectively accounting for 95.4% of the model's predictive power.H. forskaolii thrives in the midland to highland agro-ecological zones of the region, which are characterized by low seasonal mean temperatures, especially during the wettest and driest quarters.Moreover, habitats rich in evergreen trees and mosaic herbaceous species, and loam, clay loam and sandy clay soil types are particularly conducive to the growth of H. Forskaolii.Future climate change scenarios indicate a general trend of decreasing suitable habitats for H. forskaolii.Predictions under the ssp585 scenario reveal particularly concerning outcomes, with a notable decline from 4.26% in the 2030s to 19.09% in the 2090s.This highlights the detrimental impact of climate change on bee forage availability and the overall well-being of honey bees.Moreover, the assessment of current and projected habitat suitability for H. forskaolii, provides valuable scientific insights for formulating adaptation and mitigation strategies, thereby enhancing the development of climate-resilient bee forages.
Therefore, research should focus on identifying and evaluating technologies to promote climate change adaptation and mitigation efforts for the H. forskaolii species.This includes developing resilient cultivars by investigating the physiological and genetic mechanisms that enable the species to adapt to changing climate conditions.Additionally, exploring management strategies to support H. forskaolii's adaptation is crucial for effective climate change adaptation and mitigation efforts.

Figure 1 .
Figure 1.Maps showing the occurrence records of H. forskaolii in Tigray, Ethiopia.Figures were created using ArcMap (version 10.4).

Figure 2 .
Figure 2. Matured H. forskaolii at the time of flowering (right panel) and when honeybees collecting nectar (left panel) from H. forskaolii in natural habitat (source 21 ).

Figure 4 .
Figure 4. ROC (Receiver Operating Characteristic) curves and AUC (Area Under the Curve) averages for MaxEnt with 10 repeated runs.

Figure 5 .
Figure 5.The contribution of environmental factors as determined by the Jackknife analysis of regularized training gain in MaxEnt models for H. forskaolii.The variables mentioned, namely aez, bio17, bio4, bio8, and bio9 correspond to the Agroecological zone, Precipitation of the Driest Quarter, Temperature Seasonality (Co-efficient of Variation), Mean Temperature of Wettest Quarter and Mean Temperature of the Driest Quarter, respectively.

Figure 6 .
Figure 6.The relationship of major environmental factors with habitat suitability for H. forskaolii.The variables mentioned, namely bio4, bio8, bio9, and aez correspond to Temperature Seasonality (Co-efficient of Variation), Mean Temperature of the Wettest Quarter, Mean Temperature of the Driest Quarter and Agroecological zone, respectively.The areas under the response curve for agro-ecology, agroecological zones (AEZ 1 to 6), and land uses (10 to 120) are detailed in the legends of Figs.7c and 8a.

Figure 7 .
Figure 7. Honey bee population density (a), annual precipitation (b) and Agroecological zonation (c) and mean annual temperature (d) in Tigray.To determine the colony population density, we have used average data (2004-2021) from the Central Statistical Agency of Ethiopia (CSA).The honey bee population density is expressed as the average number of honey bee colony population per square kilometer.Figures were created using ArcMap (version 10.4).

Figure 8 .
Figure 8. Land use types (a) and Soil texture (b) that influence the distribution of H. forskaolii in Tigray.The numbers in the legend, ranging from 10 to 120 in the land use classification (i.e., based on the Land Cover CCI (Geoportail UCL-Geomatics) and from 1 to 9 for soil texture types, respectively used for modeling the species distribution and produce their response curves.Figures were created using ArcMap (version 10.4).

Figure 9 .
Figure 9. Current habitat distribution of H. forskaolii (left panel) and its suitability status (right panel) in Tigray.Figures were created using ArcMap (version 10.4).

Table 1 .
List of variables considered in predicting the distribution of H. forskaolii: bioclimatic variables (temperature and precipitation), topography, land use and soil conditions.

Table 3 .
Relative contribution of environmental variables against existing climatic conditions as explained by the Maxent model.